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Abstract 


The scattering length for the 77 -meson collision with deuteron is calculated 
on the basis of rigorous few-body equations (AGS) for various r/N input. The 
results obtained strongly support the existence of a resonance or quasi-bound 
state close to the rjd threshold. 

PACS numbers: 25.80.-e, 21.45.+V, 25.10.+S 

The production of 77 mesons and their collisions with nuclei have been studied experimen¬ 
tally and theoretically with increasing interest during the last years. To a large extent this is 
motivated by the fundamental questions of charge-symmetry breaking and the break-down 
of the Okubo-Zweig-Iizuka rule. Another relevant question concerns the possible formation 
of 77 -nucleus quasi-bound states. 

In many respects the 7 /-meson is similar to the 7 T°-meson despite it’s being four times 
heavier. Both are neutral and spinless, they have almost the same lifetime, rv./ 10 18 sec, 
and are the only mesons which have a high probability of pure radiative decay, that is, 
their quarks can annihilate into on-shell photons. However, when being involved in nuclear 
reactions they behave rather differently. The S’n-resonance (1535), for instance, is formed 
in both 7 tN and rjN systems, but at different collision energies, 

E™fi(Su) = 1535 MeV — tun — m n ~ 458 MeV 
E^n(Sh) = 1535 MeV — m N — m v ~ 49 MeV . 

Thus, due to the large mass of the 77 -meson (547.45MeV), this resonance is very close to 
the rpY-threshold. Furthermore it is very broad, with T se 150 MeV, covering the whole 
low energy rjN region. As a result the interaction of nucleons with 77 -mesons in this region, 
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where the S'-wave interaction dominates, is much stronger than with pions. 

After its creation the A r *(1535)-resonance decays into ijN and nN channels with equally 
high probabilities [l|] 


( N + rj (35 - 55 %) 

iV*(1535) -f \ N + n (35 - 55 %) (1) 

[ other decays (< 10 %) , 

which indicates that the rjN and nN interactions are to be treated by a coupled channel 
analysis. The resulting rjN interaction, obtained in this way, turned out to be attractive 
[2J. This raises the question whether the attraction is strong enough to support ^-nucleus 
quasi-bound states. Let us recall in this context that, because of their short lifetime, 77 - 
mesons can only be observed in final states of certain nuclear reactions. Within nuclei they 
are considered to undergo multiple absorption and production processes via the Su reso¬ 
nance, with a final transition into pions. Such quasi-bound states therefore would be of 
considerable interest for studying 77 -meson properties in more detail. 


For the calculation of these states various model treatments were employed, among them 
the optical potential method |3|-|(|, the Green’s function method [0 , and the modified multi¬ 
ple scattering theory [|8|. Calculations, based on the exact Alt-Grassberger-Sandhas (AGS) 
equations [H, for the 7 ]d system were also made in Ref. 


10 in early nineties. 


The predictions concerning the possibility of 77 -mesic nucleus formation are very diverse. 
One obvious reason for such a diversity is the poor knowledge of the r]N forces. Another 
reason comes from the differences among the employed approximations some of which might 
be detrimental in view of the resonant character of the tjN dynamics and the delicacy of the 
quasi-bound state problem. As was shown in Ref. 0 this problem cannot be adequately 
addressed by a meson-nucleus optical model or any low-order perturbation theory. 


Among the approximate approaches the few-body dynamics of the 77 -nucleus systems 
was most explicitly treated in our previous calculations [ |TT| -|H>|| based on the finite-rank 
approximation (FRA) of the nuclear Hamiltonians. The shortcoming of these calculations 
is the neglect of excitations of the nuclear ground states. This appears justified in the 77 4 He 
and possibly in the 77 -triton ( 3 He) case, but is quite questionable in 77 -deuteron collisions. 


I 11 the present paper we, therefore, treat the 77 -deuteron system on the basis of the exact 
few-body equations (AGS). Both the NN and t]N amplitudes entering them are chosen 
in separable form which reduces the dimension of these equations to one. The same r/N 
amplitude have been used in the FRA calculations. This allows us to compare our present 
calculations of the 7 ]d scattering length with the previous approximate results, i.e., to ex¬ 
amine the effect of the neglect of nuclear excitations employed in the FRA. It turns out 
that the discrepancies are not large for most of the i]N parameter sets. This indicates that 
the conclusions drawn in our previous investigation [|TT| -jT5f were already fairly reliable and 
should be even more reliable in the less sensitive 77 -triton or 77 4 He cases. 
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The 7/-deuteron scattering length is the value of the elastic scattering amplitude 

/(Pi, Pi; *) = — (2 tt) 2 /ii (pi; I t/n (^) I pi; (2) 

at zero collision energy. Here the subscript 1 labels the rj(NN) partition and the 77 -deuteron 
channel whose asymptotic states are normalized as 

(pi;V’d|Pi;V’d> = <Kp'i - Pi) • 

The transition operator Uu obeys the system of AGS equations 

3 

Uij(z) = (1 - 8 ij )gQ l {z) + - 5 ik )t k (z)g 0 (z)U kj (z) , i,j = 1,2,3 , (3) 

k =1 

where go is the free Green’s function in the three-body space, and t t the two-body T-matrix 
for the i-th pair (ti = t^N)- For both t vN and t NN we used one-term separable forms 

ti(z) = \Xi)n(z){Xi\ ■ (4) 

For the NN subsystem Eq. (dj) implies that the asymptotic wave function is related to the 
form-factor yq) according to 


IPu'fpd) = ^o(^)IXi)IPi) , ( 5 ) 

at z = p\j 2/i! + E d with E d being the deuteron energy. Due to 0) and the scattering 
amplitude (||) can be rewritten as 

/(Pi, Pi! 2 ) = -(2vr) 2 /ii(p'il^n(-)|Pi) , (6) 

where the operators X lv defined as 

x ij(z) = (Xi\9o(z)U ij (z)g 0 (z)\x j ) , 


obey the system of equations 


Xij(z) = Zij(z) + £ Z ik (z)r k ( z - X kj (z) 


k =1 


2 g k 


(7) 


with 


Zij(z) = (1 - Sij)(Xi\go(z)\xj) ■ 


The identity of the nucleons implies that X 31 = X 2 i, r 3 = r 2 , and Z 3i = Z 2 \ , which reduces 
the system to two coupled equations 


p\ 


Xn(z) = 2 Z 21 (z)r 2 I * - XX \ x 21 {z) , 


v\ 


X 2 i(z) — Z 2 i(z) + Z 21 (z)t 1 I z — J X n (z) + Z 23 (z)t 2 I z 


vl N 

2/^2 


( 8 ) 


X 21 (z) . 
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Eventually, after making the S-wave projection of the matrix elements (p'lRRjpj) and 
(p'|Zjjjpj), we end up with one-dimensional integral equations which can be solved numer¬ 
ically by replacing the integrals by Gaussian sums. 


The S'-wave separable nucleon-nucleon and 77 -nucleon T-matrices of the form (0 were 
adopted from Refs. |16| and JTTJ]. ffowever the parameters originally proposed in Ref. |T6 
for the T-matrix 


t N N{p',P',z) = — n(p')- 

v(p) = WTW ’ 

A(z ) = — tgh ( 1 




A{ Z )B{z_ 


-v(p) , 


r 00 

B(z) = / — 

Jo z — 


z 

E c 

„ 2 „, 2 , 


p 2 v 2 (p) 


p 2 /m N + is 


dp , 


were slightly modified to E c = 0.816 fm -1 , /3 = 1.604 fm -1 , and y 2 = 1.883 fm -2 which 
correspond to more recent values of the triplet NN scattering length, o/vjv = 5.424 fm, and 
the effective range r^N = 1.759 fm ||17| , |I8| . With these parameters the deuteron is bound at 
2.205 MeV and has an R.M.S.-radius \/< r 2 >d = 1.887 fm. 


Instead of treating r/N and nN as a two-channel system it is customary to describe 
the t/N interaction by a one-channel complex potential. The strength parameter A of the 
corresponding T-matrix, 


t V N(p'iP-,z ) = 


A 


Or / 2 + a 2 )(z -E 0 + iT/2){p 2 + a 2 ) ’ 
is chosen to reproduce the complex 77 -nucleon scattering length a V N, 

a\E 0 - zT/2) 


(9) 


A = 


(27r) 2 /i^7v 


~drqN 


( 10 ) 


the imaginary part of which accounts for the flux losses into the nN channel. The range 
parameter a in (||) is fixed in a somewhat more complicated way (see Refs. [2j.|20|l ), while E$ 
and T are the parameters of the SR resonance |l|], 

E 0 = 1535 MeV — (mjv + m v ) , T = 150 MeV . 

The two-body scattering length a V N, which defines the strength parameter A via (|10| ) , is 
not accurately known. Different analyses |19|] provided for a V N the values in the range 


0.27 fm < Re a r/i v < 0.98 fm , 0.19 fm < Irn a V N < 0.37 fm 


( 11 ) 


The parameter a is also known with large uncertainty. Three different values are given in 
the literature, namely, a = 2.357 fm -1 p, a — 3.316 fm -1 p0fl , and a = 7.617 fm -1 0. We, 
therefore, calculate the 77 -deuteron scattering length A vd for values of a^jy and a covering 
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these intervals. The results of our calculations are given in Tables 2, and 3, and also shown 
in Figs. 1 and 2. 


In order to check our numerical procedure, we perform test calculation of the scattering 
length with decreasing values of the meson mass, and compare the results obtained with the 
corresponding scattering lengths given by 


a fsa _ 

A t)d ~ 


d 


7rcr 


r j 

d 

Jo 

it a 4 r 


e~ ar 3 

r + - (3 + ar) -- 


i(£ , o^r)j Mr)| 2 dr, 


where Ud is the radial wave function of the deuteron. This formula is easily derived in the 
Fixed Scatterer Approximation (FSA). As it should be, the AGS and FSA results converge 
to each other when the target particles become much heavier than the incident one (see 
Table 1.) . 


In Table 2 we compare the present AGS calculations with our previous results obtained 


by means of the FRA [O]. In Table 3 we present A v d calculated with a = 3.316 fm 1 for 


various values of a V N given in the literature. For comparison, we show also the results of 
three different approximate calculations: that of Ref. || where two versions of the Multiple 
Scattering Theory (MST) were used, and a new FRA-calculations which we performed with 
the deuteron wave function (^|). In contrast to our previous FRA calculations this wave 
function (which in the coordinate representation is of the Hulthen form) provides the same 
NN input as in the AGS calculations. The dependences of A vd on Rea^ and Ima^jv are 
shown in Fig. 1 and Fig. 2. 


The curves depicted in Fig. 1 are similar to Argand plots, though they represent the 
scattering amplitude as a function of the coupling constant instead of the collision energy. 
Despite this the circular movement of the points on these curves is of the same nature as 
in the genuine Argand plot. Indeed, to draw an Argand plot one moves the point on the 
energy axis from left to right in the vicinity of a resonance pole. It is clear that if we fix 
the energy instead and move the resonance pole itself from right to left, the behavior of 
the amplitude should be similar to the Argand circle. An increase of Re a V N makes the pN 
interaction more attractive which moves the pd resonance poles towards negative energies, 
i. e. from right to left. The Argand-likc shape of the curves in Fig. 1 implies therefore that 
at a certain value of Rea^ (within the interval from 0.25 fm to 1 fm) the resonance pole 
bypasses (from below) the point E = 0 and becomes a quasi-bound pole. 


It should be emphasized here that, in contrast to the genuine Argand plot, all the points 
depicted in Fig. 1 correspond to the same energy, E = 0, and therefore the range in which 
the //-deuteron S'-matrix pole moves on the energy plane, when Re a v n varies within its un¬ 
certainty interval, cannot be inferred from these circular curves. They, however, definitely 
indicate that such a pole exists and crosses the threshold line Re A = 0. The positions of 
this pole for different values of Re a V N were explicitly calculated in a previous publication 
I5| , where the FRA-approximation was employed. Recent measurements of p production 
in the reaction p + n —> d + 77 show a substantial enhancement of the cross section near 
threshold, as compared to what is expected from phase space analysis [|26|, implying the 
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existence of such a pole. 

As can be seen in Table 3, both MST and FRA fail to give the correct (especially its 
real part) in the case of strong r]N interaction (when Rea^ > 0.5 fm) while for small values 
of Rea^Ar these methods work reasonably well. Their failure in the case of strong two-body 
forces might be due to poor convergence of the multiple scattering series and to increased 
influence of the break-up channel. 


To summarize, in the present work we perform exact AGS calculations for the t]d scat¬ 
tering length for various r]N input that include new data which appeared since the first 
calculation in 1991 Jl0j. The results obtained with these new data strongly suggest that 


a resonance or quasi-bound state could exist near the r]d threshold, in agreement with the 
prediction of Ref. 
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TABLES 


TABLE I. Convergence of the AGS and FSA results for decreasing sequence of the meson 
mass values. The parameters of the rjN- potential are fixed by a V N = (0.75 + *0.27) frn and 
a = 2.357 fm" 1 . 


77 -mass 

Affi s (fm) 

( fm ) 

m^ 

3.941 + *6.702 

1.936+ *3.162 

m v / 2 

1.548 + *0.596 

1.374 + *0.856 

m v / 3 

0.891 + *0.283 

0.878 + *0.439 

m.,,/4 

0.629 + 70.185 

0.640 + *0.292 

m v / 5 

0.487 + *0.138 

0.503 + *0.218 

m v /10 

0.230 + *0.061 

0.242 + *0.095 

m^/20 

0.113 + *0.029 

0.119 + *0.045 

O 

CO 

0.075 + *0.019 

0.079 + *0.029 

m^/40 

0.056 + *0.014 

0.059 + *0.022 

7n v /50 

0.045 + *0.011 

0.047 + *0.017 


TABLE II. Comparison of i]d scattering lengths (in fm), obtained using the AGS and FRA 


methods, for 9 combinations of the parameters of the rjN- potential. 



a = 2.357 (fm" 1 ) 

a = 3.316 (fm x ) 

a = 7.617 (fm x ) 

0"qN (fm) 

AGS 

FRA 

0.71 + *0.79 

0.66 + *0.82 

0.71 + *0.84 
0.65 + *0.85 

0.71 + *0.92 

0.62 + *0.89 

0.27+ *0.22 

AGS 

FRA 

0.79 + *0.68 

0.75 + *0.73 

0.81 + *0.73 
0.74 + *0.76 

0.83 + *0.81 

0.72 + *0.81 

0.28+ *0.19 

AGS 

FRA 

1.81 + *2.44 

1.53+ *2.00 

1.64 + *2.99 

1.38+ *2.15 

0.75 + *4.00 
1.14+ *2.22 

0.55 + *0.30 






































TABLE III. Results of the AGS and three different approximate calculations of with 
a = 3.316 fm -1 . 


^IV-input 

exact A 7ld (fm) 

approximate A^ (fm) 

Ref. 

a vN (fm) 

AGS 

MST I |] 

MST II @] 

FRA 

© 

0.25 + *0.16 

0.73 + i0.56 

0.66 + i0.71 

0.66 + i0.58 

0.65 + i0.70 

n 

0.27 + i0.22 

0.71 + i0.84 

0.57 + i0.97 

0.64 + *0.81 

0.59 + i0.96 

m 

0.291 + i0.360 

0.38+ *1.36 

0.17+ *1.35 

0.42 Til.25 

0.21 + *1.35 

i 

0.30 + i0.30 

0.61 +*1.22 

0.39+ *1.28 

0.58 Til.11 

0.42 + *1.27 

0 

0.430 + i0.394 

0.50 + i2.07 

0.14+ *1.91 

0.65+ *1.73 

0.24 + *1.88 

i 

0.44 + i0.30 

1.15+ *1.89 

0.63+ *1.93 

1.01 +*1.50 

0.68+ *1.86 

m 

0.46 + i0.29 

1.31 +*1.99 

0.72 + i2.04 

1.11 Til.54 

0.76 + *1.96 

0 

0.476 + i0.279 

1.49+ *2.06 

0.81 + i2.15 

1.22 Til.56 

0.84 + i2.05 

0 

0.51 +*0.21 

2.37+ *1.77 

1.48+ *2.31 

1.65 + *1.39 

1.38+ *2.22 

i 

0.55 + i0.30 

1.64+ *2.99 

0.61 + i2.73 

1.40 + *1.98 

0.69 + *2.51 

0 

0.579 + i0.399 

0.34 + i3.31 

—0.13 + i2.64 

0.93 + i2.41 

0.13+ *2.52 

0 

0.62 + i0.30 

1.80+ *4.30 

0.36 + i3.36 

1.65+ *2.41 

0.55 + *2.95 

0 

0.876 + i0.274 

-8.81 Ti4.30 

-2.76 + i4.24 

2.42 + i5.55 

—0.67 + i3.98 

0 

0.888 + i0.274 

—8.63 + i3.49 

—2.90 + *4.12 

2.37 + 75.79 

—0.73 + i3.99 

0 

0.98 + i0.37 

—4.69 + *1.59 

-2.75 Ti2.77 

—0.06 + i6.20 

-1.18+ *3.59 
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FIGURES 



FIG. 1. The values of A v d calculated for Ima^jv = 0.30 fm while Reoy/v is changing from 
0.25 fm to 1 fm with the step 0.01 fm. An increase of Rea^Tv moves the points in the anti-clockwise 
direction along the curve trajectories which correspond to three choices of the range parameter a. 
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